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We use electromagnetically-induced transparency (EIT) to probe the narrow electron-spin reso- 
nance of nitrogen-vacancy centers in diamond. Working with a multi-pass diamond chip at tem- 
peratures 6-30 K, the zero-phonon absorption line (637 nm) exhibits an optical depth of 6 and 
inhomogenous linewidth of ~30 GHz fuU-width-at-half-maximum (FWHM). Simultaneous optical 
excitation at two frequencies separated by the ground-state zero-field splitting (2.88 GHz), reveals 
EIT resonances with a contrast exceeding 6% and FWHM down to 0.4 MHz. The resonances 
provide an all-optical probe of external electric and magnetic fields with a projected photon-shot- 
noise-limited sensitivity of 0.2 V/cm/\/Hz and 0.1 nT/\/Hz, respectively. Operation of a prototype 
diamond-EIT magnetometer measures a noise floor of <1 nT/v^Hz for frequencies above 10 Hz and 
Allan deviation of 1.3±1.1 nT for 100 s intervals. The results demonstrate the potential of diamond- 
EIT devices for applications ranging from quantum-optical memory to precision measurement and 
tests of fundamental physics. 



Electromagnetically-induced transparency (EIT) is 
an optical coherence effect which provides exquisite con- 
trol over the absorption and dispersion in atomic media. 
In atoms with two coherent ground-state levels that can 
be optically excited to the same excited state (a "A sys- 
tem"), EIT results in ultra-narrow transmissive spec- 
tral features, with resonance quality factors exceeding 
10^^ [1-4] and more than 10'' reduction in optical group 
velocity [2, 5]. Numerous EIT-based applications are 
being pursued, including precision measurement [6-9], 
few-photon nonlinear optics [10-13], optical buffers [14- 
16], and quantum optical memories [17-19]. 

Critical to EIT-based applications is the simultane- 
ous presence of substantial atomic absorption and long- 
lived ground-state coherence [20]. This has motivated 
the use of atomic gases, where high optical depth (^1) 
and long coherence times (3>1 ms) can be simultane- 
ously realized [17, 19]. However, a solid-state approach 
is desirable for compatibility with large-scale fabrication 
processes. EIT in various rare-earth doped crystals has 
been observed [21-23], and there has been substantial 
progress towards applications in quantum information. 
However, one drawback of many rare-earth-doped sys- 
tems is weak optical transitions [17, 18]. 

Ensembles of nitrogen- vacancy (NV) centers in dia- 
mond may provide an ideal compromise, owing to the 
relatively strong NV-light coupling [24, 25], which can 
be enhanced using optical microcavities [26-29], and 
long ensemble spin coherence time (3>100 ms using 
decoupling techniques [30]) at temperatures r<100 K 
[31, 32]. EIT in diamond was observed before [33], 
but there the use of a high-defect-density diamond and 
the necessity of a large magnetic field (^0.1 T) limited 
the range of possible applications. More recent studies 
showed that A systems can be realized near zero mag- 
netic field [34-38]. 

In this Letter, we report EIT with low-defect-density 
diamond at zero magnetic field. Using electron- 



irradiation and annealing techniques [39] to enhance the 
NV absorption coefficient, a, and a multi-pass diamond 
chip to increase the optical path length, L, we realize 
optical depth, aL > 1, and coherence time, > 1 /is. 
EIT provides a means of probing the NV ground-state 
spin resonances without microwave irradiation, and we 
show that high-sensitivity, all-optical electric and mag- 
netic field sensing is possible even in the presence of 
large bias electric fields. These are desirable features 
for cryogenic applications such as fundamental physics 
experiments [40, 41] and studies of novel superconduc- 
tors [42] . Our results are well-described by a model for 
EIT in inhomogenously-broadened media and can be 
extended to applications in nonlinear optics and quan- 
tum information. 

The apparatus is depicted in Fig. 1(a). Light 
(637 nm) from an external-cavity diode laser (ECDL) 
was passed through an electro-optic phase modulator 
(EOM), producing sidebands with a tunable sideband- 
carrier detuning near 2.88 GHz. The EOM output was 
combined with a green repump beam (532 nm), neces- 
sary to reverse optical bleaching [33, 34, 43]. The light 
beam then traveled 8 times through a multi-pass di- 
amond chip housed in a continuous-flow liquid-helium 
cryostat. The transmitted light was spectrally filtered 
and detected by a photodiode. Experiments were per- 
formed at TwlO K. 

The diamond chip was a chemical-vapor-deposition- 
grown, single crystal with dimensions 4.5x4.5x0.5 mm^ 
and nitrogen density [N]<1 ppm. The sample was ir- 
radiated with 2 MeV electrons (dose: 4 x 10"^^ cm~^) 
and subsequently annealed at 800° C for several hours. 
This resulted in [NV"]=25±15 ppb, as measured by op- 
tical spectroscopy [39]. Two opposing (100) sides of 
the diamond were polished and coated with Ag to pro- 
duce 80-nm-thick mirrors. On one side, two windows 
(~0.5 x 0.5 mm^) were left uncoated to permit optical 
transmission. 
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FIG. 1. (a) Experimental apparatus. Microwave frequency 
modulation (FM) and lock-in detection were used for mag- 
netometry (Fig. 4). AOM-acousto-optic modulator. The 
dark square on the diamond surface is from a (disconnected) 
gold electrode, (b) Transmission spectrum of a weak opti- 
cal probe, (c) NV level structure and allowed A transitions 
at _B = in the moderate-strain regime. '^^N quadrupole 
splitting is omitted. 



This configuration allowed 8 passes through the di- 
amond, limited by the angular deviation between the 
polished sides. An optical micrograph of the laser- 
induced fluorescence is shown within Fig. 1(a). The 
beam ('-^lOO fim diameter) was coUimated over the 
entire propagation length. Figure 1(b) shows the 
zero-phonon line transmission spectrum of a 500 nW 
probe. The optical depth reaches aL—6 with fuU-width- 
at-half-maximum (FWHM) -30 GHz. Throughout, 
we normalize transmission by its off- resonant, room- 
temperature value, to account for interface losses, and 
the laser detuning, A^, is relative to absorption maxi- 
mum (470.480 THz). 

Figure 1(c) illustrates the NV level structure. The 
center possesses Csy symmetry, with a paramagnetic 
{S=l) ground state and two S=l excited-state or- 
bitals [44, 45]. Under moderate transverse strain, 
(5_L=15±10 GHz, level anticrossings in the lower excited- 
state orbital (Er^) mix electron spin projection, permit- 
ting optical transitions from both ground-state nis=0 
and ?7is=±l manifolds [36, 38]. The TOs=0 and 
■ms~±l manifolds are split at zero magnetic field by 
-D=2.88 GHz, and hyperfine coupling with the ^^N nu- 
cleus (/=1) results in three separate A schemes. 

Figure 2(a) shows the timing of optical pulses used to 
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FIG. 2. (a) Timing diagram of optical pulses used to probe 
EIT resonances, (b) Transmission spectrum as the sideband- 
carrier detuning was swept through two-photon resonance, 
(c) EIT spectra at different magnetic fields applied along a 
[100] direction. The spectra are offset for clarity. Overlayed 
are Lorentzian fits (see text). In (b,c) Al~ — 15 GHz. 



observe EIT. Alternating green and red pulses was nec- 
essary to simultaneously achieve high EIT contrast and 
minimize bleaching. To efficiently excite both arms of 
the A transitions [Fig. 1(c)], the sinusoidal optical phase 
modulation was set to yield a sideband:carrier:sideband 
intensity ratio of 0.7:1:0.7 (Supplementary Information, 
SI). Higher-order sidebands contributed < 8% of the 
total intensity and are neglected in our analysis. 

Figure 2(b) shows the transmitted red light as a func- 
tion of EOM drive frequency, i^eom- The broad anti- 
hole is due to optical pumping [34]. When Jj^eom — i^j ^ 
^'^nati where Ai/'nat^lS MHz [46] is the homogenous 
excited-state linewidth, NV centers resonant with one 
of the excitation frequencies can be excited from one 
spin sublevel, but are eventually trapped in the other 
sublevel, resulting in high transmission. As veom ap- 
proaches 2.88 GHz, one of the sidebands acts as a re- 
pump for the carrier, so there are no trap states, re- 
sulting in lower transmission. The FWHM of this fea- 
ture (typically 50-200 MHz) depends on several factors, 
including excited-state dephasing [31], optical power 
(Prod), and spectral diffusion [27, 38], but it is always 

As t'EOM matches exactly the ground-state split- 
ting (two-photon resonance), the transmission in- 
creases sharply. These resonances are the hallmark 
of EIT, and they exhibit much narrower FWHM, 
A:/oit=0.4-1.3 MHz < Ai^nat- They occur because NV 
centers are optically pumped into a "dark" coherent su- 
perposition of ground-state levels, \D), which cannot in- 
teract with the light due to quantum interference. This 
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can be understood by considering a simplified model 
for the NV center consisting of two ground-state lev- 
els, ms=l (|1)) and ms—O (|0)) driven optically to a 
single excited state, {E^). The Hamiltonian under the 
rotating wave approximation is: 

Jf=h/{4Tr){no\E,){0\+ni\E,){l\) 

+ h{A,-Ao)\l){l\ + hA,\E.,){E,\ + h.c, ^ ' 

where h is Planck's constant, and ilg and Ag are, re- 
spectively, the Rabi frequency and detuning of the \s) -O- 
\Ex) transition. On two-photon resonance (Ai=A2), 
the state, |_D) = ^===(f2i|0)— f^oll)), is completely 

decoupled from the optical fields, satisfying J^\D)=0. 
The orthogonal superposition is coupled, so NV centers 
are pumped into \D), resulting in increased transmis- 
sion. 

The EIT linewidth is limited by the decoherence rate 
of the dark superposition, 7^^, and its narrow width al- 
lows for sensitive, all-optical probing of the NV ground- 
state level structure. Figure 2(c) shows EIT spectra 
for several values (B) of magnetic field along a [100] 
direction. [100]-directed fields preserve the degeneracy 
of the four NV axes, enabling higher EIT contrast. At 
B—0, the outermost resonances are split by ^2Ahf, 
where Ahf= — 2.17 MHz is the longitudinal hyperfine 
coupling constant [47]. The small splitting (5o~0.5 MHz 
between the two innermost resonances (771^=0) arises 
from transverse crystal strain, which behaves as an 
ensemble-averaged effective electric field [48]. This 
effective electric field is \E±\=Sf)/{2dgs^)f^l5 kV/cm, 
where dgs^—l-7 Hz/V/cm is the ground-state transverse 
dipole moment [49]. 

The values of veom on two-photon resonance are ap- 
proximately equal to the ground-state transition fre- 
quencies [50]: 



u,^=D±yJ{gfiBB cose + m.AHFY + {dgs^E^^, (2) 

where (7=2.003 is the electron-spin g-factor [51], 
/iB=13.996 GHz/T is the Bohr magneton, and 
cos0«l/\/3 is the projection of the field along each NV 
axis. Equation (2) neglects AC Stark shifts of order 
f^g/(167r2Ae)<10 kHz, where Ae«0.5 GHz is a typical 
frequency separation between excited states (SI). 

We fit the spectra to a sum of six Lorentzian pro- 
files (plus an offset C), with amplitude A, and FWHM 
Afcit. The central frequencies were constrained by 
Eq. (2); D and E± were fit globally, while B, Avcit, 
A, and C were allowed to vary between spectra. Reso- 
nances that satisfied {gfisB cos9 + miAHF\ < dgs^\E±\ 
(innermost resonances within the top three spectra in 
Fig. 2) were given separate amplitude (Aq) and FWHM 
(Albeit, 0) from those where the condition was not held 

(Ai, Al^eit.l) [50]. 

To gain further insight into the EIT lineshapes, we 
determined the contrast and FWHM as a function of 
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FIG. 3. (a,b) Power dependence of zero-field EIT FWHM 
and contrast (A_l~ — 13 GHz) determined from Lorentzian 
fits as in Fig. 2(c). (c,d) Photon-shot-noise-limited magnetic 
and electric field sensitivity [Eq. (3)] , inferred from zero-field 
EIT spectra, as a function of laser detuning. 



Prcd at B=0, Fig. 3(a,b). For mi=0 resonances, the 
zero-field contrast is here defined as i?o =^0/(^0 + C), 
and for rrii—ztl resonances it is Ri—2Ai/{Ai + C). 
These data were fit by a 3-level density-matrix model 
for EIT in inhomogenously-broadened media (SI) 
[34]. The fit parameters include a Rabi frequency 
conversion factor, Psat=7rAf^jjjPi.cd/^^o=2.4±l.l mW, 
the ratio ni/f^o=0.08±0.02, which reflects the 
ensemble-averaged asymmetry in A transition 
strengths [38], and nuclear-spin-dependent ground- 
state dephasing [50] 7gs.„j.=±i/(27r)=240±92 kHz and 
7gs,mi=o/(27r)=99±30 kHz. The apparent saturation 
of the contrast arises from photo-ionization (SI), which 
reduces [NV] at a rate cxP^^^ [43]. Above T^IO K, the 
contrast falls off sharply with temperature, becoming 
neghgible at T«30 K. This effect (SI) is well described 
by a nine-level model of the NV center that includes 
temperature-dependent excited-state dephasing [31]. 

Near BsaO, the EIT resonances provide a means 
to simultaneously sense electric and magnetic fields. 
The theoretical sensitivity of an electric or magnetic 
field sensor based on an optically-detected signal, S, 
is given by the minimum detectable field which gives 
signal :noise=l; ^E^nin^ ^^gf^^^^ and 5Brain= \j^sfdB\ ' 
where 5S is the standard deviation of S [52]. From 
Eq. (2), we see that dS/dE^O except when \giiBB -\- 
raiAHF\'^dgs^^\E±]. At BsaO, this condition is satis- 
fied for the rrii—Q resonances, whereas high-sensitivity 
magnetometry is possible using the mi=±l resonances 
provided \E_L\<\AHF\/dgs^fvl30 kV/cm. 

If SS is limited by photon shot noise, the magnetic 
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sensitivity is (SI): 



gfiB cos Ri y Ftm 

where P is the detected optical power, Ep is the pho- 
ton energy, and tm is the measurement time. A similar 
expression can be found for sensitivity to electric fields 
within diamond (SI), substituting dgg^Ro for gfiBRi 
[50] (neglecting crystal-strain inhomogeneity) . These 
equations assume detected fields lie along a [100] di- 
rection, but arbitrarily-oriented fields can be detected 
using a suitable bias field [53, 54]. 

We studied the projected electrometer and magne- 
tometer performance by analyzing zero-field EIT spec- 
tra as a function of A/,, Fig. 3(c-d). In both cases, 
the sensitivity is optimized in the range Ai= — 5 to 
—20 GHz. Not surprisingly, this corresponds to the 
range of \Ex) strain shifts where A transitions have been 
observed in single-NV experiments [38]. The optimal 
sensitivities (Ai= — 13 GHz) are 105±6 pT/a/Hz and 
0.24±0.04 V/cm/v^Hz for magnetometry and electrom- 
etry, respectively. Throughout we use the metrology 
convention that 1 Hz measurement bandwidth corre- 
sponds to tm = 0.5 s. 

High-sensitivity operation of our device as an all- 
optical magnetometer was accomplished using lock- 
in detection of a resonance peak [55] in a bias 
field B=Q.2 mT [Fig. 4(a)]. A microwave signal 
i'EOM=2.8717 GHz was frequency-modulated at 89 kHz 
[56], with deviation 300 kHzpp, and sent to the EOM. 
The photodetector signal was demodulated using a lock- 
in amplifier, and the in-phase signal was sent to an oscil- 
loscope. Figure 4(b) shows the magnetometer response 
when a separately-calibrated test field was applied. Ad- 
ditional calibrations are presented in SI. Summarizing, 
the magnetometer response remains linear over a range 
of > 1 ^T, and the bandwidth covers ~ 100 Hz (limited 
here by electronic filtering). 

To characterize the sensitivity, the magnetic noise 
was measured without a test field. Figure 4(c) shows 
the noise-equivalent magnetic field spectrum. The noise 
floor is <1 nT/-\/Hz for frequencies above 10 Hz. The 
floor was dominated by lock-in amplifier input noise, 
as evidenced by its persistence when the photodetec- 
tor was unplugged. The expected shot noise, based on 
Eq. (3) after incorporating the finite quantum efficiency 
of the detector, is 0.3 nT/v^Hz. 

Near 1 Hz, the noise fioor rises. A possible cause is in- 
stabilities in Picd lead to fiuctuations in AC Stark shifts 
(SI). Nevertheless, the magnetometer recovers high sen- 
sitivity for long integration times, as evidenced by the 
Allan deviation [47] plotted in Fig. 4(d). For a gate 
time of 100 s, the Allan deviation is 1.3±1.1 nT. 

In comparing with existing technologies, we are not 
aware of another sensor that simultaneously measures 
low-frequency electric and magnetic fields with high 
sensitivity. Miniature vapor-cell magnetometers [57] 
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FIG. 4. (a) Schematic of the frequency modulation (FM) 
technique employed for magnetometry. (b) Lock-in signal 
when a 16-Hz oscillating field (350 nTpp) was applied along 
a [100] direction, (c) Magnetic noise spectrum of the EIT 
magnetometer. The spike at 60 Hz is a real magnetic signal 
arising due to operation in an unshielded environment (SI). 
Also shown is the noise spectrum when the photodetector 
was unplugged, (d) Allan deviation of a 500 s data set. 

achieve low-frequency sensitivity of ^0.01 pT/\/Hz at 
T«470 K, but the sensitivity rapidly degrades with de- 
creasing temperature [52]. At low temperature, minia- 
ture superconducting quantum interference devices 
have excellent sensitivity, 5i?min^0.1 pT/-\/Hz. How- 
ever they typically suffer from 1/f noise, so DC mea- 
surements require external calibration [58]. In addition 
to its long-term stability and dual electric/magnetic 
field sensitivity, our sensor is probed entirely optically, 
in large electric fields, and it does not produce fields of 
its own. 

Future applications may benefit from several im- 
provements. Higher EIT contrast is possible using spec- 
tral filtering to detect only one optical frequency, tai- 
loring the absorption spectrum using holeburning tech- 
niques [59], or employing other magnetic/electric field 
geometries [33, 38] . The coherence time can be extended 
by orders of magnitude using dynamic decoupling [30] 
and/or with isotopically-pure diamond [47, 60, 61]. Fi- 
nally, integration with optical cavities [26, 28, 29] will 
allow higher optical depth and larger optical intensities, 
while decreasing the device volume. 

In summary, we observed narrow EIT resonances in 
a multi-pass diamond chip. The high optical depth and 
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narrow inhomogenous linewidth enable new diamond- 
based applications including all-optical electrometry 
and magnetometry. We operated a prototype diamond- 
EIT magnetometer with sub-nT/-\/Hz sensitivity and 
excellent long-term stability. Integration with photonic 
networks may enable new applications, including few- 
photon nonlinear optics and quantum-optical memories. 
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Supplementary Information: Electromagnetically-induced transparency in a diamond 
spin ensemble enables all-optical electromagnetic field sensing 
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FIG. SI. Contrast of zero-field mi=±l EIT resonances as 
a function of sidebandxarrier intensity ratio. 
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FIG. S2. Peak transmission of zero-field mi—iLl EIT reso- 
nances as a function of Pred along with quadratic fit. 



SIDEBAND POWER OPTIMIZATION 

In order to determine the optimal sidebandxarrier 
intensity ratio, we varied the microwave power driv- 
ing our EOM and fit the B=0 EIT spectra to obtain 
the contrast. Figure SI shows the EIT contrast versus 
the first-order sidebandxarrier intensity ratio, at con- 
stant Prod=l-3 mW. The sidebandxarrier ratios were 
determined using a scanning Fabry-Perot interferom- 
eter. For ratios below 1, the higher-order sidebands 
contribute ^10% to the total intensity and can be ne- 
glected. For ratios above 1, the higher-order sidebands 
become substantial and lead to a reduction in contrast, 
as the intensity for any given frequency is lower. In 
the experiments presented in the main text, we chose 
a sidebandxarrier:sideband intensity ratio of 0.7:1:0.7, 
where the maximum contrast was observed. 



MODELING EIT RESONANCES AT T<10 K: 
3-LEVEL MODEL 

At temperature T<10 K, the excited-state levels of 
the NV center are well resolved [SI], with typical fre- 
quency spacing, Ae«0.5 GHz. Consequently, for rela- 
tively weak optical excitation, fio,! << SttAe, we can 
treat our system as an ensemble of three-level atoms. 
The situation at higher temperature, where phonon in- 
teractions lead to significant broadening of the excited 
states, is treated in the next section. 

We begin by writing the rotating-wave Hamiltonian, 



as in Eq. (1) in the main text, but now in matrix form: 

/ na/2 \ 
Jf/h = 27r(5 rJi/2 , (SI) 
\ no/2 ni/2 27rAo / 

where the basis is {|0), {Er^)}. Here is the de- 
tuning from the |s) O l-E^) resonance ("one-photon de- 
tuning"), (5=Ao — Ai is the detuning from two-photon 
resonance, and we have assumed real Rabi frequencies 
ri*=17s. Note that at zero field, the ms=±l resonances 
are degenerate, so in principle the system involves four 
levels. However only one superposition of ms—±l can 
be optically coupled to \Ex), so we label this superposi- 
tion |1) and ignore the uncoupled, orthogonal superpo- 
sition. Unlike the EIT dark state, defined in the main 
text, this uncoupled superposition does not play any 
role in the EIT spectrum, as it does not include |0). 

In order to compute steady-state solutions, we use a 
master equation for the density matrix, p, in the pres- 
ence of relaxation, ^[p]: 

dp/dt = ~i[.Jf/h, p] + M[p] = 0. (S2) 

Here ^ includes spontaneous emission, Az^nat=15 MHz, 
as well as nuclear-spin-dependent ground-state decoher- 
ence, 7gs,mi, and longitudinal ground-state relaxation, 
7i. We set all other dephasing terms to zero. 

We further make the assumption that the medium is 
optically thin. This may seem like an incorrect assump- 
tion, as aL ranges between 1.3-3 (Fig. S2). However we 
believe the assumption is justified because most of the 
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absorption is incoherent background absorption from 
NV centers which do not have A transitions and/or be- 
long to hyperfine states which are far detuned from two- 
photon resonance. Under this approximation, the EIT 
contrast and width can be calculated by evaluating the 
\Ex) population, P\e^),\e^)j as a function of two-photon 
detuning, using Eq. (S2). In order to incorporate the ef- 
fect of inhomogenous broadening of \Ex) levels, we sum 
over spectra from variable Aq over a range —1 to 1 GHz 
(much larger than fig) [S2]. 

In fitting the linewidth and contrast data [Fig. 3(a,b)] 
to spectra calculated in this way, we make the following 
assumptions, fls is constrained to be oc \/ Prcd ■ Igs.mi 
is allowed to be different for resonances associated with 
different nuclear spin projections (TOi=0 and mi=±l). 
This is justified because the mi=0 and nii—ztl hyper- 
fine levels are sensitive to different noise sources (elec- 
tric versus magnetic field noise, respectively), as de- 
scribed in the main text and in [S3]. 71 is constrained 
by the observed antihole width [Fig. 2(c)], W (in Hz), as 
7i=A:y„atf^g/(27rVF2) [S4]. The value of 71/ (27r) is then 
of order a few kHz. This may seem to contradict recent 
results which predict 7i=l/Ti«27r x 0.01 Hz [S5]. How- 
ever, those results were for NV centers in the absence 
of optical fields. We find that the in the presence of op- 
tical fields, the effects of, for example, spectral diffusion 
and interaction with other excited-state levels lead to a 
much larger effective 71 that describes our results. 

The EIT contrast is sensitive to incoherent absorp- 
tion from NV centers which do not exhibit A transi- 
tions, hyperfine states which are off two-photon reso- 
nance, and the removal of NV~ centers due to photo- 
ionization. The photo-ionization rate, in particular 
changes the power dependence, as it is quadratic in 
Pred [S6]. Figure S2 shows the peak transmission of 
EIT spectra as a function of Pred- The quadratic fit 
is further evidence that photo-ionization plays an im- 
portant role in the EIT spectra. To account for these 
effects when fitting the simulated contrast to experi- 
mental data, we modify the calculated EIT contrast, 
Pcaic, as P„,=am,Pcaic/(l + bP^^^), where am, and b 
are fitted parameters. 

Under these approximations, we fit the model 
to the data in Fig. 3(a,b) and find the following 
fit parameters: Psat=7i'Ai'^a^4Prod/f^o=2.4±l.l mW, 
17i/rJo=0.08±0.02, 7g^,,„^=±i/(27r)=240±92 kHz, 
7ss,m,=o/(27r)=99±30 ' kHz, a„,=±i=1.0±0.1, 
a„,=o=0.3±0.1, and 6=0.4±0.2 mW^^ 



MODELING EIT RESONANCES FOR HIGHER 
TEMPERATURE: 9-LEVEL MODEL 

As briefiy mentioned in the text, the EIT contrast 
falls off rapidly with increasing temperature. In our 
experiments, temperature is measured via a diode lo- 
cated on the copper "cold finger" on which the dia- 
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FIG. S3. (a) Temperature-dependence of the 

experimentally-determined contrast for the mi=±l zero- 
field EIT resonance (r2o=27r x 7.8±2.8 MHz). Also shown 
is the simulated contrast, using the 9-level model, for three 
values of fio- The simulated contrast is reduced by a factor 
of 13 to match the low-temperature experimental contrast, 
(b) Experimental contrast, shifted by 7 K, along with cor- 
responding simulations. The simulated contrast is reduced 
by a factor of 13(35) for TOi=±l(0) resonances. 



mond is mounted. Figure S3(a) shows the EIT contrast 
of mi=±l resonances at B=0. The contrast becomes 
negligible for cold-finger temperatures T>30 K. In gen- 
eral, this phenomenon is due to phonon-induced orbital 
relaxation the excited state [SI]. 

For the highest power used here, Piod=l-3 mW, the 
Rabi frequency can be calculated, using the 3-level- 
model fits described above, as f2o=27r x 7.8±2.8 MHz. 
As this Rabi frequency is much smaller than the en- 
ergy spacing between levels in the excited state, the 
3-level model can in principle be used to qualitatively 
describe the observed temperature dependence. This 
can be done by including temperature-dependent de- 
phasing, TesiT) for both \s) O transitions, where 
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Tes{T) is determined using parameters in Ref. [SI]. In 
this model, the contrast is approximately proportional 
to r2§/[res(T)7gs^mJ. This would suggest that the dele- 
terious effect on T^g, from increasing temperature, can 
be overcome by increasing Qq. However, due to interac- 
tions with other excited-state levels, we find this is not 
the case. 

A more general treatment considers interaction with 
all six excited state levels. For this 9-level model 
(all three ground-state sublevels and the six excited- 
state levels), we work in the Zeeman basis and include 
excited-state spin-orbit, spin-spin, and Stark-effect in- 
teractions [S7], pure ground-state spin dephasing and 
longitudinal relaxation (determined using 3-level model 
fits above), and excited-state homogenous broadening, 
Afnat=15 MHz. We also make the following assump- 
tions: 

• inhomogenous broadening is neglected. 

• the transverse strain is 15 GHz (near a level ani- 
crossing) . 

• all optical coupling frequencies, Qq, are equal. 

• the laser frequencies are resonant with the lowest 
excited state. 

• phonon-induced excited-state relaxation affects 
only the orbital portion of the wavefunction and 
preserves m^. 

• the contrast is calculated using the excited-state 
population (optically-thin medium). 

Figure S3(a) shows the calculated contrast for three 
different values of Qq. The simulated contrast is re- 
duced by a factor of 13 to match the low-temperature 
experimental contrast. We find that increasing i7o to 
27r X 100 MHz only increases the contrast by a small 
amount, extending the operating temperature range by 
^ 10 degrees. At higher values of Hq the FIT lineshapes 
are severely distorted due to interaction with other ex- 
cited states. Moreover, for f7o > Stt x 100 MHz, no 
significant improvement in contrast at any temperature 
is observed. Note that in our experiments, the maxi- 
mum Rabi frequency was riosa27r x 8 MHz. At these 
low excitation rates, both the 3-level and 9-level models 
predict that the contrast is still increasing with increas- 
ing power, as observed in Fig. 3(b). 

We also find that better quantitative agreement be- 
tween experiment and theory can be obtained if we as- 
sume that the local temperature in the beam path is 7 K 
greater than what is measured at the cold finger. Figure 
S3(b) shows the shifted experimental values along with 
their corresponding theoretical curves. The change in 
local temperature may be due to a combination of sub- 
stantial optical absorption (^10 mW, time-averaged) 
and relatively poor thermal contact. 




FIG. S4. Magnetic response, 13(B), as a function of bias 
field along a [100] direction. 



PHOTON-SHOT-NOISE LIMITED 
SENSITIVITY AS A FUNCTION OF B 

Equation (3) in the main text describes 
the photon-shot-noise-limited minimum de- 
tectable field, 5Bmim under the assumption that 
\dS/dB\—g^BCosORi//S.Vcit- Here the contrast is 
Ri=2Ai/{Ai + C) for 5=0, but it is about half, 
A/{A + C), when a sufficiently large bias field lifts the 
ras—±.l degeneracy. We define the higher-field slope as 

More generally, the slope, dS/dB, depends on both 
the lineshape and the bias field. It can be expressed as: 

\dS/dB\ = qP{B)^. (S3) 

The factor of g is a numerical constant of order unity 
that accounts for the lineshape; q— 3 4=1.290 for a 
Lorentzian lineshape. The bias-field-dependent factor, 
(3{B), arises due to the overlap of Ams=±l resonances 
near B=0. This factor is important because it can dic- 
tate which bias field is optimal for a given value of Ai'eit. 

Figure S4 shows a plot of P{B), assuming Lorentzian 
lineshapes. At exactly zero magnetic field, /3(0)=0 
and the conventional lineshift magnetometry em- 
ployed here is not possible. As B increases, f3{B) 
rises sharply, reaching a maximum /3{Bopt)=2 at 
g^slBopt cos6'|=At'eit/(2\/3). For large fields, satisfy- 
ing gfiB\B\^Ai'eit/2, P{B) approaches unity. 

The sharp dependence of dS/dB, and consequently 
SB^im on bias field may have important implica- 
tions for proposed fundamental physics experiments. 
For example, Ref. [S8] proposes to search for the 
electric dipole moment of neutrons in an ambient 
magnetic field B=l ^T. Using typical experimen- 
tal parameters in this work (P=100 ^W, A/{A + 
C)=0.035, Ai/cit=l MHz, cos6l=l/^/3, (7=3^3/4), we 
find that, at |i?|=l /iT, the minimum detectable field 
would be 6B„iin—i20 pT/y/Hz. This is nearly an or- 
der of magnitude larger than the optimal sensitivity, 

dBrmn-50 pT/^ at |B„pt| = 18 /iT. 
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However, improvement in sensitivity at low bias 
field ((//j,B|i?|<cAi'oit/2) can be realized by decreasing 
Ai'oit- This is because in this regime /3(i?)ocl/Ai/cit, 
so consequently i5SmmOcAz/^;^ [Eq. (S3)]. If instead 
Ai'cit=0.25 MHz, and the contrast does not change, 
then SB„,in^27 pT/VSz at ^T, a 16-fold im- 

provement. Narrowing of the EIT resonances could be 
accomplished by reducing Pj-cd or using isotopically pure 
diamond. Optimal sensitivity could also be realized by 
actively applying an appropriate bias field. 

We note that, as discussed in the text, electric-field 
sensitivity is realized only under sufficiently low bias 
magnetic field. For rrii—O resonances, the sensitiv- 
ity falls off by a factor of ~-\/2 under a bias field, 
\B„iax\=dgs^\E±\/{gfiB cos 9). For the effective field 
|i?j^|~15 kV/cm present in our sample, this value is 
|^mo2:|~15 /iT. Fortunately, this bias field is com- 
parable to the bias field where SBj^m is optimal, 
I -Bopt 1=18 AtT, ensuring that high-sensitivity electrome- 
try and magnetometry can be simultaneously realized. 



COMMENT ON E-FIELD SENSING 

The electric-field sensitivity quoted here is for fields 
inside of diamond. If the goal of the electrometer is 
to sense external charges/fields, one needs to account 
for dielectric screening. In most geometries, the electric 
field measured in vacuum will be larger than the field 
in diamond by a factor of order the dielectric constant, 
5.7. 



MAGNETOMETER LINEARITY 

In the present implementation, the range over which 
the magnetometer output remains linear, ASiin, is con- 
strained by the EIT linewidth. In particular, using 
the microwave frequency-modulation technique, ABun 
is always smaller than the modulation depth, /S.Vawd- 
Ai^niod is chosen to optimize the magnetometer scale 
factor (conversion of lock-in output to magnetic field 
units), Ai^mod ~ Ai/oitZ-v/S (see Fig. S4 and related 
discussion). 

Figure S5 shows the lock-in signal as a function of car- 
rier frequency using the experimentally-optimized mod- 
ulation depth, Afniod = 300 kHzpp. From the fig- 
ure, we see that the scale factor remains linear in a 
range of At^un « 100 kHz, corresponding to A£?iin = 
Ai^iin/(gMB cos 61) w 20 

To verify the linear response, we applied 16-Hz sinu- 
soidal current modulations to a coil with different ampli- 
tudes, and measured the resulting EIT magnetometer 
output. The conversion of applied current to applied 
-Bioo-field was extrapolated by measuring the shifts in 
EIT spectra to large changes in current (as in Fig. 2(c) 
in the main text). Figure S6(a) shows the lock-in re- 
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FIG. S5. (a) Lock-in signal as a function of microwave 
carrier frequency. 
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FIG. S6. (a) EIT magnetometer response to 16-Hz ap- 
plied fields of three different amplitudes, (b) Measured field 
amplitude as a function of the expected applied field ampli- 
tude. The solid line corresponds to perfect linearity (slope 
= 1) with no offset. 



spouse to each applied field. The measured amplitude 
is compared to the separately-calibrated input field am- 
plitude in Fig. S6(b). As expected, the magnetometer 
remains linear for the two smallest applied fields, cor- 
responding to a range ASun ~ 1 /xT. For the largest 
amplitude (1.1 /iTrms) there was some deviation from 
the expected response; this may be due to, for example, 
a mismatch of the lock-in reference phase. 



MAGNETOMETER SCALE-FACTOR DRIFT 

The long-term accuracy of the sensor relies on the 
scale factor remaining constant over time. To test the 
scale-factor drift, we applied a 16-Hz sinusoidal field for 
50 s, and monitored the measured amplitude over time 
by fitting a sine curve to 0.5-s intervals. Figure S7 shows 
the results of the measurement. No clear long-term 
drift of the measured amplitude is observed, so we con- 
servatively bound the scale- factor drift at < 1%/min. 
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FIG. S7. (a) Fitted amplitude of a 16-Hz applied field over 
the course of time. 
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FIG. S8. (a) Attenuation of the EIT magnetometer signal 
due to the electronic filtering used here. 

Note that any residual drift can be corrected for by, for 
example, periodically monitoring the lock-in response 
curve or adding an additional EOM frequency modula- 
tion with a different depth and frequency. 

MAGNETOMETER BANDWIDTH 

The fundamental limit on our EIT magnetometer 
bandwidth is set by the polarization rate of NV cen- 
ters into the dark state, ~ Afoit. In the present imple- 
mentation, the bandwidth was limited by the electronic 
filtering of the lock- in amplifier. The amplifier used a se- 
ries of four 6 dB/octave low-pass filters, resulting in an 
attenuation to signals at frequency, /, of (l-|-27rrL/)~^. 
In order to avoid artifacts introduced by the 4 kHz rep- 
etition rate of the optical pulses, a filter time constant, 
= 1 ms, was used. 

Figure S8 shows the electronic attenuation as a 
function of frequency. We tested the magnetometer 
frequency response by varying the frequency of a ^ 
130 uTj-ms sinusoidal field modulation, and found ex- 
cellent agreement with the response plotted in Fig. S8. 
This effect is also responsible for the slight downward 



slope of the noise floor at frequencies > 10 Hz in Fig. 
4(c) of the main text, as those data were not corrected 
for the small frequency-dependent attenuation. 

MAGNETOMETER NOISE 

From Figures 4(c-d) in the main text, it is evident 
that there is low-frequency noise near 1 Hz. We investi- 
gated the origin of this noise by recording noise spectra 
under different conditions, including: 

1. the photodetector was unplugged, and we 
recorded the lock-in output [Fig. 4(c)]. 

2. t'EOM was detuned —3 MHz from two-photon res- 
onance and we recorded the lock-in signal with 
the photodetector connected. 

3. the signal from a commercial fluxgate magnetome- 
ter, placed close to the position of the diamond, 
was recorded. 

However, in all these cases, the noise was absent. It is 
worth mentioning that the fluxgate magnetometer also 
picked up magnetic noise at 60 Hz of a similar amplitude 
to the noise peak recorded by the EIT magnetometer 
[Fig. 4(c) in the main text]. However the noise floor 
near 1 Hz was <1 uT/a/Hz. 

We estimate that AC Stark shifts of order 
riQ/(167r^Ae)~10 kHz, are present under our typical 
operating conditions. This was confirmed by EIT spec- 
tra obtained from the nine-level model described above. 
Consequently fluctuations in Picd of 1% lead to varia- 
tions in magnetometer output of order a few nT. Such 
fluctuations are consistent with the magnitude of the 
noise. We note that fluctuations in power also cause 
the overall transmission to change, but this generally 
happens slow compared to our modulation frequency 
(89 kHz) so it does not show up in the magnetome- 
ter signal (as evidenced by the lack of the noise when 
t'EOM was detuned —3 MHz). The AC-Stark fluctu- 
ations may be caused by vibrations due to turbulent 
cryostat flow and/or vibrations of the transfer line con- 
necting the cryostat to the liquid helium dewar. Fu- 
ture low-frequency applications may benefit from active 
power stabilization or post-processing the magnetome- 
ter signal based on monitoring the power. 
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